Four Fermion Operator Matching with NRQCD Heavy and AsqTad Light Quarks 
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We present one-loop matching coefficients between continuum and lattice QCD for the heavy-light 
four-fermion operators relevant for neutral B meson mixing both within and beyond the Standard 
Model. For the lattice theory we use nonrelativistic QCD (NRQCD) to describe b quarks and 
improved staggered fermions (AsqTad) for light quarks. The gauge action is the tree-level Symanzik 
improved gauge action. Matching to full QCD is carried out through order as, A.qcd / Mh, and 
asl{aMb). 
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I. INTRODUCTION 

— mixing, with q = d or q = s, has been the fo- 
cus of much attention by experimentalists and theorists 
in recent years, and it will continue to be probed exten- 
sively in the LHC era [l| . On the experimental front the 
mass differences between the "heavy" and "light" eigen- 
states in the B° systems, AM^ and AMg, are now known 
very accurately Q . Measurements of the decay width 
difference AFg and the phase 4>s by D0 and CDF have 
also appeared [1, H, @]. On the theory side neutral B 
systems are of particular interest as a possible window 

into New Physics (NP). In the Standard Model -B° - B° 
mixing does not occur at tree level and must go through 
box diagrams involving the exchange of two FF's at low- 
est order. NP could enter through the exchange of new 
particles in the box diagrams, or through new tree level 
contributions. Studies of the neutral B meson param- 
eters can impose important constraints on different NP 
scenarios 

Theoretical estimates of mixing rates employ effective 
Hamiltonians involving four fermion operators. Matrix 

elements of these operators between B^ and B, states are 
needed to complete the calculations and this requires con- 
trol over non-perturbative QCD. For instance the Stan- 
dard Model expression for the mass difference AMq is 
given by , 



AM, = ^j^\ 



V;qVtb\'v2Soixt)MBj'BBB,, (1) 



where Xt = m^/Myy, 77^ is a perturbative QCD cor- 
rection factor, So{xt) the Inami-Lim function and Vtq 
and Vtb the appropriate Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. The nonperturbative QCD in- 
put into this formula is the combination /^^ 13 b^ where 

/b is the Bq meson decay constant and Bb the renor- 
malization group invariant bag parameter. Lattice QCD 
provides a first principles approach to obtaining these 
crucial non-perturbative factors. The first realistic lat- 
tice results that include effects from two very light and 



the strange sea quarks (so-called Nf = 2 -f- 1 simulations) 
have now appeared ^, Ji, Jil] . 

An important step in all lattice calculations is the 
matching between four-fermion operators in continuum 
QCD that enter into formulas such as eq.([T]) and the op- 
erators of the lattice theory used in the lattice QCD sim- 
ulations. In this article we present one loop matching 
results for a complete basis of AB — 2 four fermion op- 
erators relevant for mixing both within and beyond the 
Standard Model. In the lattice theory we employ a non- 
relativistic QCD (NRQCD) action for the b quarks, an 
improved staggered quark action (AsqTad) for the s and 
d quarks and the Symanzik improved glue action. 

A subset of these results was used already in reference 
Q in the calculation of all the hadronic matrix elements 
relevant for the determination of AMs and AFg in the 
Standard Model. That was the first lattice determination 
of these parameters with Nf = 2-1-1 sea quarks. 



II. THE FOUR FERMION OPERATORS AND 
MATRIX ELEMENTS IN QCD 

In order to study neutral B meson mixing phenom- 
ena, we focus on the following five AB — 2 four-fermion 
operators (i and j are color indices) 
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The subscript q stands for either the d or the s quark, 
both of which we take to be massless in our matching 
calculations, and Pr^l = i 75). 

Operators Ql, Q2 and 03 appear in the Standard 
Model and are relevant for the mass and width differ- 
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FIG. 1: Diagrams depicting the two Dirac structures SI = 
{uq'^ i,2Uq)(vQY' 2, iVq) aud S2 = (uQTi^2Vq){vQV2,iUq). The 
incoming and outgoing states are those of eq.([9}. "A", "B", 
"C", and "D" are color indices. 




ences AM^ and AF^. Matrix elements of Ql, for in- 
stance, lead to /^^ Bb^ ■ Additional AB — 2 operators, 
such as Q4: and Q5 are required when going to extensions 
of the Standard Model. The above five operators go un- 
der the name of the "SUSY basis of operators" in the 
literature [HI, . At intermediate stages of the match- 
ing calculation, we found it useful to introduce matrix 
elements of two more operators. 



Q6 = [^h''PL%) {^il.PR% 

Q7 = (W'.-rPL^) {^1.Pr% 



(7) 
(8) 



Within perturbation theory matching can be carried 
out by considering scattering from an incoming state with 
[heavy antiquark + light quark] into an outgoing state 
with [heavy quark + light antiquark] . Symbolically, 



(9) 



where "A", "B", "C", and "L>" are color indices. We 
also introduce external Dirac spinors Uq and Vq for the 
incoming light quark and outgoing light antiquark respec- 
tively, and similarly uq and vq for the outgoing heavy 
quark and incoming heavy antiquark. 

At the tree-level, matrix elements of operators Ql, Q2, 
QA, and Q6 become, 

{0Ut\ (^iri*;) [%r2^^q) \in)tree = (10) 
+ 5ab5cD [ {uQTlUq){vQV2Vq) + {uQV2Uq){vQVlVq) ] 
- 5ad5cB [{uQTlVq){vQT2Uq) + {uQV2Vq){vQVlUq)] . 

Diagrammatically, we will denote the two Dirac struc- 
tures , SI = {uQTi^2Uq){vQT2,iVq) aud 5*2 = 
{uQTi,2Vq){vQT2^iUq) by the first and second diagrams 
in Fig.l respectively. 

For operators QS, Qb and Ql one has similarly, 



{out\ {^iT.^l) {^iT2%) \in)t 



(11) 



FIG. 2; One-loop corrections to four fermion operators Qk. 
The four external lines are the same as in Fig.l. 



- 5ab5cD [ {uQVlVq){vQT2Uq) + {uQT2Vq){vQT lUq) ] 
+ 5adScB [ {uQTlUq){vQT2Vq) + (jlQT2Uq){vQVlVq) ] . 

For the cases where Fi = r2 the two terms in the square 
brackets in egns. pO] ) and (|lip are identical and can be 
combined with a factor of 2. When Fi 7^ F2 these terms 
should be kept separate during the matching procedure. 
The RHS of these equations can often be simplified using 
Fierz relations. In particular, using formulas of Appendix 
A one finds. 
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(12) 



and hence QQ and Q7 are not independent operators. 
Nevertheless, we have found the bookkeeping to be sim- 
plified if one first projects onto {QQ)tree and {Ql)tree and 
then uses cq. p^ to relate back to QS and (^4. 

Continuum one-loop corrections to tree-level matrix el- 
ements are obtained by evaluating the diagrams in Fig. 2a 
- Fig.2f'. We carry out these calculations in the M5-NDR 
scheme with definitions of general dimension four-fermion 
operators (i.e. of "evanescent" operators) as given in Ap- 
pendix B [H, [3 . We retain only terms of 0{as) and dis- 
card 0{ash.QCD/M) contributions, where M is the mass 
of the b quark. We thus expand the full QCD calculations 
to the same order to which our effective theory (lattice 
NRQCD) perturbative calculations of the next section 
are carried out. In this limit the heavy quark spinors 
obey, 



uq7o = UQ, 



(13) 



As is well known, there is mixing among the four-fermion 
operators at one-loop. For instance. 



(Ql)^^ = {Ql)tree 

-a, "cn(Ql)il + Ci2(g2)(°L 



(14) 



3 



All matrix elements are taken between {out\ and \in) 
and the superscript (0) means that we are working with 

spinors obeying (fT3|) . Similarly, one finds that {Q2)^''^ 
has contributions from Q2 and Ql, {Q3)^^^ from Q3 and 
Ql, (g4)^ from Q4 and Q6, and (Q5)^ from QS and 
QT. The one-loop coefficients Cxy depend on the b quark 
mass M, the gluon mass A that acts as an IR regulator, 
and the AfS'-NDR renormalization scale fj,. Their values 
for the operators in the basis ([2]) are given by, 



antiquarks {^-q)- At lowest order in 1/M, matrix ele- 
ments such as (fTO)) and pT|) . are reproduced by working 
in the effective theory with 



If one introduces an effective theory field, 
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— e// — 

then and the QCD field ^'b are related by a Foldy- 
Wouthuysen-Tani transformation. In particular. 



(29) 



C22 



C21 



1 

47r 



16, 4 ^2 



1 [3 1, 2 A2 

4^ 2 + 3^"m^ + 3^"aP 



(17) 
(18) 



where the V acts to the left. By inserting (|29p into the 
expressions for the four-fermion operators Ql - Q7 , one 
sees that 0{Aqcd/M) corrections to (l?7)) can be ob- 
tained within the effective theory by adding the following 
1/M operator corrections 
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eq. (|T2|) C46 and C57 can be replaced by. 



The one-loop coefficients cn, C12, C22 and C2 1 ag ree with 
the values already published in reference Using 



(25) 
(26) 
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III. MATRIX ELEMENTS IN THE EFFECTIVE 
THEORY 

In effective theories such as HQET or NRQCD one 
works separately with heavy quark fields that create 
heavy quarks (^'q) and with those that annihilate heavy 



Ojl 



1 

2M 



v^Q ■ 7ri*,) {^Qr2^g) (30) 
(v% • 7r2*g)] + (Ti r2). 



As mentioned already, our effective theory consists of lat- 
tice NRQCD for heavy quarks, the improved staggered 
AsqTad for light quark action and Symanzik improved 
glue action. Details about these actions and their Feyn- 
man rules can be found for instance in reference [16\ . In 
this effective theory, we have evaluated one-loop correc- 
tions to matrix elements of O and Ojl. The one- loop 
corrections to O involve the same diagrams Fig. 2a - 2f 
as in continuum QCD. One obtains 



{Qiyff = [l + a,cf,]{Ql)f,l + a.cf2(Q2)i°L, 



(31) 



and similarly for all the other Qfc's. Several of the lat- 
tice one-loop integrals are IR divergent and we use a 
gluon mass A to extract the IR finite contributions as 
explained in reference [l^. The c^y's have the same IR 
divergent ln(A2) terms as the corresponding Cxy in con- 
tinuum QCD, and they also depend on the bare heavy 
quark mass. The divergent terms will cancel when we do 
the matching (see below) and things are reduced to finite 
differences such as [cxy — c^y]- We carry out the lattice 
perturbative calculations in both Feynman and Landau 
gauges and use gauge invariance as a check on our results. 

In order to calculate the one-loop renormalization coeffi- 
cients for the matrix elements of the operators Ojl the 
diagrams of Figs. 3 & 4 need to be evaluated. One finds. 



{Qijir" - (giji>il 



has 



11 



:ree 

(32) 




FIG. 3: One loop corrections to the 1/M corrections Ojl. 
The dark dot denotes a derivative acting on either the heavy 
quark or the heavy anti-quark propagator. We only show 
corrections associated with Fig.2a and 2b. Similar diagrams 
exist for corrections to 2c - 2f'. 



where we ignore 0{as) corrections to {Q\jl)f^^^. Similar 
expressions are obtained for the other Ojl. 

The coefhcients C,^^ tell us about the "mixing down" 
of a dimension seven operator Ojl onto dimension six 
operators O. On dimensional grounds these coefficients 
go as l/(aM), "a" being the lattice spacing. They rep- 
resent "power law" contributions from matrix elements 
in the effective theory. Power law terms are unavoid- 
able when working with effective theories and they need 
to be subtracted in order that the effective theory pro- 
duce the same physics as full QCD (which does not suffer 
from power law contributions) . Our matching procedure, 
described in the next section, will be such that power 
law contributions are removed from matrix elements of 
Ojl through 0{as/{aM)). Errors from the mismatch 
between QCD and effective theory due to power law con- 
tributions will come in at 0{a1/{aM)). 



IV. 



MATCHING 



We wish to relate the continuum QCD matrix elements 
(Qk)'^'^^ to the matrix elements {QkY^^ in the effective 
theory. The latter will ultimately be replaced by out- 
put from nonperturbative simulations. We will focus on 
matching of Ql. The other Qfc's are handled identically. 

If one expands the first term on the RHS of in 

powers of 1/M, i.e. {Ql)tree (Ql) 
then this equation becomes 



(0) 



(0) 



MS 



[1 + as cii](Ql)il 
as ci2(Q2)|°^, 



(33) 



{Qmirl- 



The next step is to rewrite the matrix elements {■■■)tree 
appearing on the RHS in terms of matrix elements in 
the effective theory {...Y^^. This can be accomplished 
by inverting (pij) and ([5^ . To the order that we are 




FIG. 4: Further possible diagrams associated with 1/M cor- 
rections Ojl. In practice, however, these do not contribute 
to the (^^''. Diagrams where the derivative acts on the heavy 
anti-quark propagator (top right line) similarly do not con- 
tribute. 



working one has, 

(Qi)il = {Qiyff 



.{Q2)\ 



(0) 

ree 



,(Q2) 



(34) 
(35) 



and 

(Qiji 



,(0) - 



{Qiji 



(36) 

Upon inserting the last three equations into (|33p one ends 
up with, 



(Ql)"" = 

[1 + as pn] {Qiy" 



(37) 



+ 0{alasAQCD/M), 



C'' {Q2r"' 



where, 



Pxy — C-xy C^y (38) 

With (I37p we have achieved the goal of relating the con- 
tinuum full QCD matrix element {Ql)'^^ to matrix ele- 
ments in the effective theory. Although this perturbative 
matching calculation was carried out with external scat- 
tering states, one carries over matchings such as l\'S7\i to 
the case of hadronic matrix elements between and 
B'^q states where then {QlY^f or {QljlY^f must be 
evaluated nonperturbatively. With such nonperturbative 
matrix elements in mind we define 



{Qljiyff -as [Cm"^^ + C^{Q2) 
= {Qljl 



eff- 



(39) 



' sub ■ 
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The combination (Qljl)^^^ represents thus the matrix 
element of the dimension seven 1 /M correction in the ef- 
fective theory with power law contributions subtracted 
out through 0{as/{aM)). This matrix element gives 

us the physical Aqcd/M contributions to (Ql)*^'^ up 

2 

to corrections of 0{ ^"^^ ). Further discussion of power 
law subtractions for the case of heavy-light currents are 
given in reference [17|. A more complete derivation of 
one-loop matching formulas including all contributions 
at 0{asAQCD/M) is provided (again for heavy-light cur- 
rents) in reference J^]. 

A final technical detail is that eq. ([37|) differs from 
eq.(lO) of [1| in that here we assume the same normal- 
ization of states in the effective theory as in continuum 
full QCD. For the purposes of evaluating matching co- 
efficients it is convenient to do so. Any differences in 
normalization of states are taken care of at the stage of 
doing the non-perturbative calculations and of extracting 
decay constants and bag parameters. 



V. RESULTS 

In this section we summarize results for the effective 
theory coefficients c^y and C"^^ of eans. (|3T|) and p2|l. and 
for the matching coefficients pxy of eg. ([57)1 . We present 
numbers for three values of the bare heavy quark mass 
in lattice units, aMp, corresponding to the b quark mass 
on lattices with spacings 0.09/to, 0.12/to and O.lT/m, 
as fixed in previous studies of the T system [l^. These 
values of a correspond to the so-called MILC fine, coarse 
and super-coarse lattices, which have been used exten- 
sively in recent studies of heavy-heavy ^19] and heavy- 
light [l^, [2l[ quantities with the same choice of lattice 
actions as the one in the present article. 

In Table I we list values for the one-loop renormaliza- 
tion coefficients c^y after subtracting the IR divergent 
In(aA)^ pieces, together with the one-loop which are 
IR finite. Table II shows values for p^y at scale fi equal to 
the heavy quark mass M. The parameter n in these two 
tables is the stability parameter of the NRQCD action. 

In Table III we illustrate how the different diagrams 
contribute to C44 for gMq = 2.8. Calculations were done 
in both Feynman = 1) and Landau = 0) gauges 
to check for gauge invariance. We have evaluated C44 
by collecting all contributions that are proportional to 
SabScd and that can be written (using Fierz relations 
to convert where necessary) in terms of Dirac struc- 
tures (uQPLUq)(vQPRVq) OT {UQ PRUq){vQ PLVq) . We 

could just as well have collected terms proportional to 
—SadScb and of Dirac structure {uQPLVq){vQPR,Uq) or 
{uQPiiVq){vQPLUq) to obtalu the same final result. Ta- 
bic IV illustrates different contributions to cfg again 
for qMq = 2.8. Here we collect terms proportional to 
SabScd and of Dirac structure {uqY PLUq)ivQ"f^PRVq) 
or {uqj'^ PRUq){vQ'-fi,P]^Vq). From cfg one easily obtains 
C45 = —2 cfg. 



aMo (n) 


1.95 (n = 2) 


2.8 (n = 2) 


4.0 (n = 2) 




-1.196 


-0.735 


-0.403 




-1.802 


-1.315 


-0.960 


C22 


0.010 


0.014 


-0.004 


C2I 


0.020 


-0.018 


-0.050 




-0.890 


-0.644 


-0.483 


i 


0.133 


0.064 


0.010 


C44 


0.692 


0.599 


0.529 


C46 


0.097 


0.029 


-0.022 




-0.886 


-0.553 


-0.311 


c| 


-0.467 


-0.383 


-0.322 


C45 


-0.194 


-0.058 


0.044 


'-•54 


0.934 


0.766 


0.644 


c" 


0.218 


0.166 


0.123 




0.874 


0.662 


0.492 




0.364 


0.276 


0.206 




0.009 


0.007 


0.005 


^33 


-0.073 


-0.055 


-0.042 




0.064 


0.048 


0.035 


^44 


0.309 


0.235 


0.175 


^46 


-0.046 


-0.034 


-0.026 




0.091 


0.069 


0.051 




0.064 


0.048 


0.035 


^45 


0.092 


0.068 


0.052 


^54 


-0.128 


-0.096 


-0.070 



TABLE I: The one-loop coefRcients c^y and for three 
values of aMo. n is the stability parameter in the NRQCD 
action. IR divergent In(aA)^ terms are omitted. Numerical 
integration errors are of order one or less in the last digit. 



VI. SUMMARY 

We have completed the one-loop matching of a com- 
plete set of AS — 2 heavy-light four fcrmion operators 
through 0{as,h.QCD/Mi,,as/{aMi,)). The main results 
arc the coefficients pxy of Table II and the Cf^^ of Table 
I. 

We find that with the lattice actions employed in this 
article (and in our simulations) matching coefficients are 
all well behaved. None of them arc particularly large, 
and in fact many are considerably smaller than one. An 
interesting feature that can be extracted from the results 
in Table III and IV, and which holds also for the ma- 
trix elements of the other operators, is that the one-loop 
matching coefficients are dominated by the current-like 
diagrams, a and h, and the wave function renormaliza- 
tions. The contributions from the pure four- fcrmion dia- 
grams are at least an order of magnitude smaller. We are 
currently investigating non-perturbative matching meth- 
ods for heavy (NRQCD)-fight (staggered) currents. The 
same methodology could be applied here for four-fermion 
operators to calculate non-perturbatively the main con- 
tribution to the renormalization coefficients. In this way 
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aMo (n) 


1.95 (n = 2) 


2.8 (n = 2) 


4.0 (n = 2) 


Pii 


0.693 


0.462 


0.357 


Pl2 


1.165 


0.678 


0.323 


p22 


0.927 


1.000 


1.094 


p2\ 


0.029 


0.028 


0.022 


P33 


0.873 


0.703 


0.618 


P31 


0.035 


0.065 


0.082 


P44 


0.628 


0.923 


1.192 


P46 


0.052 


0.077 


0.085 


P55 


0.694 


0.563 


0.520 


P57 


0.258 


0.131 


0.027 


P45 


-0.104 


-0.154 


-0.170 


P5i 


-0.516 


-0.262 


-0.054 



TABLE II: The one-loop matching coefficients pxy in (|38p 
with 11 = M for three values of aMo. 



we could considerably reduce the uncertainty associated 
with the matching process, which is one of the main 
sources of error at present in our calculation of fs, \/ Bb, 
and fs^ \fB^a- 

The matching calculation in this article is an impor- 
tant part of the HPQCD collaboration's studies of Bs and 
Bd meson mixing via lattice QCD methods both in the 
Standard Model [1, i, h2\ and beyond. Values for the 
mass and decay width differences, AMg and AFg with 

q ~ s,d, as well as for the ratio ^ = IEiV^^^ which 

fBay/BBa 

use the results presented here will be available soon [23 |. 
Extensions of this work to matching with other lattice 
actions such as the HISQ light quark action are straight- 





C = 


= 1 


C = 


= 


diagram 




^ln(aA)^x 




^ln(aA)^x 


a + b 


0.836 X 2 


-|x2 


0.3274 X 2 





a' + b' 


0.0823 X 2 





0.0823 X 2 





c + d 


-0.0495 X 2 


i x2 


0.0139 X 2 





e 


0.0331 


1 

3 


-0.0304 


1 

2 


e' 


0.0269 





0.0269 





f 


0.0640 


1 
fi 


0.0006 







-0.924 


4 
3 


-0.416 





Zq 


-0.338 


8 
3 


0.171 


-4 


Total 


0.5996 


7 
2 


0.5993 


7 

2 



TABLE III: Contributions to the coefficient C44 from the 
diagrams of Fig. 2 for aA'Io = 2.8. ^ = 1 and ^ = refer 
to Feynman or Landau gauge respectively. The second and 
fourth columns give the IR finite contributions. Columns 3 
and 5 list IR divergent terms in units of ^ln(aA)^. Zq and Zq 
are the light and heavy quark wave function renormahzations 
respectively and are taken from reference (l6|]. The last row 
gives the full c^j. 





c = 


= 1 




-- 




diagram 




4^MaA)2x 




1 




a' + b' 


-0.0743 X 2 


i x2 


0.0209 X 2 







c + d 


0.0137 X 2 





0.0137 X 2 







e 


0.0045 





0.0045 







e 


0.0497 


1 
2 


-0.0457 




3 
4 


/' 


0.0960 


1 

4 


0.0009 







Total 


0.029 


3 
4 


0.029 


3 
4 



TABLE IV: Same as Table III for the coefficient cfg. 



forward and are planned for the future. 
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APPENDIX A: FIERZ RELATIONS IN 4 D 

In this appendix we collect some Fierz relations that 
were found to be useful in our calculations. They are 
written in terms of fixed external spinors rather than 
fermionic fields, which means that interchanging two 
spinors does not bring in a minus sign. In our calcu- 
lations, fermionic signs come in at the stage of doing 
the Wick contractions, for instance to get the RHS's of 
ea. (fTU)l and pT|) . 

(Al) 

[ugPLVq] [vqPLUq] = (A2) 
- [ugPLUq] [vqPLVq] + - [uqCt'^'' PLUq] [vqa^^PLVq] 

where cr^'' = ^[7^,7'']. 

[ugPLVq] [UgPRUq] = i [uQ-f" PllUq] [vQj^PLVq] (A3) 
[uQ-f" PLVq] [vQj^PRUq] = 2 [ugPRUq] [wQ^LW,] (A4) 

In the large M limit, i.e. with spinors obeying (fT5|). the 
relation (|A2[) simplifies to [23j . 

[uqPLVq] [vqPLUq] = (A5) 

[ugPLUq] [vQPLVq] + ^ [HQj^PLUg] [TJq Pt ] ■ 
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APPENDIX B: SOME FORMULAS IN GENERAL 
DIMENSIONS 

We carry out the continuum one-loop calculations us- 
ing dimensional regularization in the MS and NDR 
scheme. As is by now well known, this information is 
insufficient to fix one's renormalization scheme unam- 
biguously. One needs to specify in addition how one 
handles d-dimensional Dirac structures appearing at in- 
termediate stages of the calculations and how they are 
projected onto some 4-dimensional basis. This is because 
the Dirac algebra is infinite dimensional for non-integer 
d, and requires in addition to the 4-d basis set an in- 
finite set of "evanescent operators". The 0{e) (we use 
d = 4 — e) terms in the projections onto the 4-d basis 
defines one's choice of evanescent operators, and is con- 
vention dependent [H, [13, [H, [2^ . In our calculations 
we have adopted the following conventions taken from 
the literature [H, ^ 

[uQ7°'l^l''PLUq] ^QlalpluPLVq] (Bl) 
^ (16 - 2e)[uQYPLUq] [vQl.PLVq] 

[ugri^YPLUg] [vQ"fu"ffno.PLVq] (B2) 

[uQl'^l^PLUq] [vQlo^lpPLVq] (B3) 

^ {^^2e)[uQPLUq\[vQPLVq\ 
-(8 - 4.e)[uQPLVq\ [vqPLUq] 



In adopting IjBSp we are following the conventions of ref- 
erence p3 ]. 

We note that [l^ has different conventions from [2^ . 
Instead of eq. (|B3P the [1^ conventions translate into the 
relation 

^Ql^'l^PLUq] ^QlalpPLVq] (B6) 
^ (8 - e)[uQPLUq] [vqPLVq] 
-8 \UQPLVq\ [ugPiUq] 

Using (jB6|) leads to different constant terms in results for 
the one-loop coefficients C22, C21, C33 and C31. The In^^ 
and InA^ terms remain unchanged, however. Denoting by 
Cij the coefficients obtained by using the ^25*1 conventions, 
one finds that eqns.([I7D, (HID, (HH) and ^ are modified 
to, 

5. = ^(6-K^ln4-iln41 (B7) 
47r \ 3 Af2 3 M2 J ^ ' 

= 4^|3 + 3^^^ + 3^^^M^| 
1 r 8 /i2 4 A2 1 



^Q^°'^^ PLUq] ^Q^aliiPRVq] (B4) Carc is required to ensure that a consistent set of con- 

^ (4 -t- 2e)[uQPiUg] ^qPRVq] ventions is applied to different parts of a calculation con- 

tributing to a physical quantity such as AM^. 

[UQ-i'^l^PLUq] [vQ^p-ic.PRVq] (B5) 
^ (4 - A€)[uQPLUq] ^qPRVq] 
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